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Grobner Bases in Perturbative Calculations 

Vladimir P.Gerdt a 

a Laboratory of Information Technologies, Joint Institute for Nuclear Research, 141980 Dubna, Russia 

In this paper we outline the most general and universal algorithmic approach to reduction of loop integrals to 
basic integrals. The approach is based on computation of Grobner bases for recurrence relations derived from 
the integration by parts method. In doing so we consider generic recurrence relations when propagators have 
arbitrary integer powers treated as symbolic variables (indices) for the relations. 
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1. Introduction 

Since its invention in |1I2| . the method of inte- 
gration by parts has become the commonly used 
platform for reduction of loop integrals to ba- 
sic (master) integrals. The distinguishing feature 
of this method is rapidly increasing volume of 
symbolic algebraic computation needed for ma- 
nipulation with recurrence relations (RR) when 
the number of internal and external lines or/and 
loops increases, especially in the presence of such 
physical parameters as masses, external momenta 
and the space-time dimension. Over the last years 
a number of algorithmic ideas, approaches and 
software packages has been reported |3I4I5I6| to 
reduce, given a set of RR, loop integrals to a min- 
imal set of master integrals. 

There are two main aspects of an algorithmic 
approach to this problem: (i) universality, i.e. ap- 
plicability to the most general form of RR and 
(ii) computational efficiency. Most of practically 
done computations are based on the approach of 
Laporta 0] since it deals directly with RR and 
provides a rather efficient reduction procedure 
based on Gaussian elimination. However, this 
elimination is restricted to RR with specified val- 
ues of indices. As specification takes place, one 
extracts only partial information from RR. More- 
over, in typical multiloop calculations it may lead 
to an enormous number of intermediate RR 6 
that can annihilate the computational efficiency. 

Another approach suggested by Tarasov |.'iloj 
appeals to potentially much more universal al- 
gorithmic techniques based on computation of 



Grobner bases (GB). But for all that, RR are 
transformed into partial differential equations. 
Such a transformation can only be done for inter- 
nal lines endowed with different masses which are 
independent variables for the differential equa- 
tions. Insertion of these extra variables may in- 
crease substantially the running time and storage 
space required for computing GB. Furthermore, 
if one considers indices of RR as (discrete) sym- 
bolic variables for the differential system, then in 
the aggregate with extra independent variables 
(masses) this tends to increase the volume of com- 
putation so much that becomes hardly possible 
for problems of interest in modern multiloop cal- 
culations. Thus, in practice, one can treat differ- 
ential systems for loop integrals only with fixed 
numerical powers of propagators. 

In this paper we consider RR generated by 
the integration by parts method as linear finite- 
difference equations for loop integrals whose coef- 
ficients can (polynomially) depend on indices and 
all relevant physical parameters such as masses, 
scalar products of external momenta and the 
space-time dimension. We present basic concepts 
of the GB technique (Sect. 2) as the most univer- 
sal algorithmic tool for RR with symbolic indices. 
GB allow one to analyze algorithmically the com- 
plete set of algebraic restrictions imposed on loop 
integrals by RR. Specifically, we apply (Sect. 3) 
the GB technique to find a minimal set of basic 
(master) integrals which are not fixed by RR, and 
to reduce any other loop integral to these basic 
ones. We illustrate this technique by a simple ex- 
ample of one-loop propagator type integral taken 
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from 3; . We also briefly discuss computational is- 
sues of the GB approach applied to RR (Sect. 4). 

The concept of GB was invented by Buchberger 
in 1965 for systems of multivariate commutative 
polynomials and then extended to systems of lin- 
ear partial differential equations, noncommuta- 
tive polynomials and other algebraic structures 
(see, for example, book Due to their al- 

gorithmic universality, GB have found numerous 
fruitful applications. Special program modules 
for their computing in the case of commutative 
polynomials built-in all modern computer alge- 
bra systems such as Maple, Mathematica and oth- 
ers. In addition, Maple has special librarian func- 
tions for computing GB for differential equations, 
which were used in |3I5| , and for noncommutative 
Ore algebras [Hj- 

2. Grobner Bases for Recurrence Relations 

Since any tensor multi-loop Feynman integral 
by shifting space-time dimension can be reduced 
to scalar integrals we shall, without the loss of 
generality, consider the following scalar integrals 

Z*= f d d k 1 ---d d k L 1 (1) 

J lli=l i 

with n internal lines. Here propagators Pi have 
symbolic exponents vi which we collect together 
into the multi index v = \y\, v 2l ■ ■ • , v n }. RR de- 
rived from integration by parts form a homoge- 
neous system of linear finite-difference equations 
in the integral Z„ as function of its indices. Let 
Vi — \i be the minimal value of the i-th index en- 
tering in the RR system. Here Xi > are explic- 
itly known integers. The set of all multi indices 
with nonnegative components will be denoted by 
N^ , and, hence, A £ N£ . 

Let /i = {v x - \x,V2 - A 2 , • • . , v n - A„} and 
the integral I,, ^ be a function of indices in 
v. Then the RR are a set of equalities fj = 
including finite sums of the form 

f 3 =Y,KD a l^ j = l,...,p (2) 

Here D a = D" 1 ■ ■ ■ D% n with a = {a x , . . . , a n } £ 
N" Q , and each Dj is the right-shift operator 
for the i-th index, i.e., DiX^ = 2 Ml ,..., Mi +i,..., Al „. 



Coefficients W a depend polynomially on in- 
dices {v\, . . . ,v n ~} and physical parameters, e.g., 
masses, scalar products of external momenta, 
space-time dimension d. Hereafter we shall call 
integrals D 01 !^ and sums fj in J5J by (difference) 
terms and (linear difference) polynomials, respec- 
tively. Coefficients V a of polynomials fj will be 
considered as elements in the field Q{v) of ratio- 
nal functions in indices whose coefficients in turn 
are parametric rational functions. 

Consider now the set of all linear polynomials 

<F>={^VD"(/)|/eF, bp e Q(y) } (3) 

generated by the polynomial set F — {fi, . . . , f p } 
defined in 

Any element g £ < F > yields the finite- 
difference equation g = which is a consequence 
of the initial RR, and < F > accumulates all the 
consequences. The set F of difference polynomi- 
als is called a basis of < F >. 

Let y be a total order on terms D a 2 ll such 
that for any a, (3, 7 £ iV> the following holds 

(i) D"!^ y 1^ -4=> a contains nonzero indices, 

(ii) D^D a T^ y D^D fi Z^ D a X^ y D^. 

Such term orders are nothing else than admis- 
sible term orders in commutative polynomial al- 
gebra ^Hli if one compares multi indices of the 
terms. As typical examples of admissible or- 
ders we indicate "lexicographical" order >~i cx with 
ct ^icx P when ot\ > Pi or when u\ = /3± and 
a 2 > Pi, etc., and "total degree" order a ^tdeg P 

when J2 a i >J2 Pi or J2 a i = J2 Pi an d a Hex P- 
Given a term order y , one can extract the lead- 
ing term lt(/) from every difference polynomial /. 
Its coefficient is called the leading coefficient and 
will be denoted by lc(/). Now we can define GB 
for J2J as a finite polynomial set G satisfying 

1. < G >=< F >. 

2. For every / £ < F > there is g £ G and 
7 £ N^o such that lt(/) = D"f(lt(g)). 

Note that conditions (i) and (ii) of admissibility 
for >- provide existence of GB and termination of 
an algorithm for its construction exactly as they 
do in commutative polynomial algebra [lUlllj . If 
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the term order >- is considered as a certain "sim- 
plicity relation" between integrals, then condition 
(i) means that integral X, is " simpler" than that 
obtained from 2^ by shifting up some indices. 
Condition (ii) shows that such a simplicity rela- 
tion between two integrals is stable under identi- 
cal shifting of their corresponding indices. 

For two multi indices a and /?, a is called a 
divisor of (3 and (3 is called a multiple of a, if 
(3 — a £ A™ . The least common multiple of a 
and (3 will be denoted by lcm(a, (3). Thus, prop- 
erty 2 means that the leading multi index of any 
consequence of RR as element in < F > has a 
divisor among multi indices of the leading terms 
in GB. Property 1 implies that both sets of RR 
defined by F and G are fully equivalent since any 
consequence of one set is a consequence of another 
set and vice versa. A term whose multi index di- 
vides (is multiple of) another term is said to be a 
divisor (multiple) of this term. 

An important notion of the GB techniques is 
reduction. Polynomial / is said to be reducible 
modulo polynomial g if / has a term u with a 
coefficient b G <9(a*) such that u is multiple of 
lt(g), i.e., u = D~i \t{g) for some 7 G N£ . In this 
case an elementary reduction step is given by 

If /* in (J3JI is still reducible modulo g, then the 
second reduction step can be performed, etc., un- 
til an irreducible polynomial obtained. By prop- 
erty (i) of term order >-, the number of elementary 
reductions is finite . Similarly, one can reduce 
a polynomial / modulo a finite polynomial set F 
by doing elementary reductions of / modulo indi- 
vidual polynomials in F. The reduction process 
is terminated with polynomial / irreducible mod- 
ulo F. In this case we say that / is in the normal 
form modulo F and write / = NF(f,F). The 
above described procedure gives an algorithm for 
computing the normal form. 

From the above definition of GB G for set F, 
it follows that NF(f, G) = for any polynomial 
/ in < F >. Therefore, if G is known, one can 
algorithmically detect the dependency of an extra 
recurrence relation r on set F . It suffices to verify 
whether NF(r, G) vanishes or not. 



It is remarkable, that there is an algorithm for 
construction of GB for any finite set of linear dif- 
ference polynomials and an admissible term or- 
der. This algorithm was discovered by Buch- 
berger in 1965 for commutative algebraic poly- 
nomials. Below we present the Buchberger algo- 
rithm in its simplest form |1UI11| adapted to lin- 
ear finite-difference polynomials. First, we define 
the S -polynomial S(f,g) for a pair of (nonzero) 
difference polynomials / and g as follows. Let 
lt(f) and lt(g) have multi indices a and /3, re- 
spectively, and let 7 = lcm(a, /?). Then 

Buchberger algorithm: 

Start with G := F. 

For a pair of polynomials /j.,/2 G G: 

Compute S(f 1 ,f 2 )- 

Compute h := NF(S{f u fo),G). 

If ft, = 0, consider the next pair. 

If h 7^ 0, add h to G and iterate. 

Inter-reduction of the output GB, that is, reduc- 
tion of every polynomial g G G modulo G \ {g}, 
gives a reduced GB. By normalizing (dividing by) 
the leading coefficients, we obtain the monic re- 
duced GB that is uniquely defined |1UI11| by the 
input polynomial set F and term order y. 

3. Generic Master Integrals 

Having computed GB for the set © generated 
by the initial set ©, one can algorithmically, just 
as in commutative algebra |llj . find the maxi- 
mal set of difference monomials (integrals) which 
are independent modulo set Q. In other words, 
there is no nonzero polynomial in @ composed 
from those monomials, and any other monomial 
(integral) can be reduced to (expressed by means 
of) these monomials. Thus, the maximal inde- 
pendent set of monomials (integrals) is exactly 
the collection of basic integrals. They may be 
called generic master integrals to emphasize their 
relevance to RR with symbolic indices. The re- 
duction of any integral to "masters" can also be 
done algorithmically (Sect. 2) by means of GB. 

By keeping track of all intermediate elementary 
reductions one can obtain an explicit expres- 
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sion of any loop integral in terms of (generic) mas- 
ter integrals and their multiples. As we defined 
in Sect. 3, a multiple of the integral is obtained by 
applying D 1 with 7 G iV> . ft is clear, however, 
that £> 7 2" M can be evaluated only if the depen- 
dence of 2 M on indices is explicitly known. 

But how to obtain the set of master integrals? 
It can be easily detected from the leading mono- 
mials of GB. Namely, all those monomials that 
are not multiple of any leading monomial in GB 
are independent (masters). This fact is very 
well known in theory of commutative Grobner 
bases and apparently applicable to linear 
finite-difference polynomials. Clearly, any mul- 
tiple of a leading term in GB is reducible to 
(expressible in terms of) a linear combination of 
monomials that are not multiple of any leading 
term in the GB and with coefficients from Q{v). 

Though the particular form of master integrals 
depends on the choice of term order y, their num- 
ber is invariant on >- and fully defined by initial 
RR. In fact, if one accepts the interpretation of 
>~ as simplicity relation (Sect. 2), the set of mas- 
ter integrals detected from GB contains the "sim- 
plest" integrals 1 . Usually a " total degree" order is 
preferable (cf. |4I6| ) from the computational point 
of view since total degree GB are typically com- 
puted much faster then lexicographical ones. 

For pure illustrative purposes consider RR for 
the one- loop propagator type integral X Vll/2 taken 
from paper [3j and put m 2 = q 2 = in those 
relations, thus, keeping space-time dimension d 
as the only parameter 

Vi^vi-i 1/2+1 — {d — v\ — 2u2)I Vl u 2 = 0, 
v \ Zui+l V3-I ~ 1^2^-1*5+1 + \V2 — v\)T VlV3 = 0. 
Their form J5J for X Vl -i tVa —i ^ is given by 
[v-lDI -{d-v x - 2v 2 )]Tu 1 -i,v 2 -i = 0, 
[viDf - v 2 D\ + (v 2 - vi)]X Vl -iu 2 -i = 0, 
where D\, D 2 are the right-shift operators for 
indices V\ and v 2 , respectively. A total degree 
Grobner basis at D\ >~ D 2 calculated by Maple 2 

1 For numerically fixed values of indices and particular in- 
tegrals sometimes one uses other simplicity relations which 
do not satisfy (i) and (ii) UJ. 

2 Algebra of shift operators is a particular case of the Ore 
algebra [8], and one can use it for computing GB for RR 
as already mentioned in However, this way of compu- 



gives the following GB form of RR 

[(d -vx- 2v 2 )D x D 2 - vx&& T Vl -xu 2 -x = 0, 
[v 1 {d-v 1 -2v 2 )Dl + 

(vi(2v 2 - v\) + v 2 {2v 2 - d))D%]l Vl _ ll/2 ^i = 0, 
D 2 2 Ul -i U2 -i = 0. 

An immediate consequence of this GB is the 
equality 2^-1 „ 2 +2 = that follows from the last 
element of GB. This implies D 1 T- vl -\v 2 +2 = 
for 7 £ Ny . Similarly, the order D 2 y D x 
gives T Vl+2u2 -i = and D^X Vl+ i V2 _ x = for 
7 G N™ that can also be verified by calculating 
the normal form iVF(D|2 i/1 _ ll/2 _i, GB) = 0. 

The leading terms in the above GB arc 
represented by the operators D\D 2 , D\ 1 D 2 
which generate four master integrals 

T-V\ — 1 1/2 — 1 1 2l/j I/ 2 — 1 ) Zj/j — 1 l/ 2 I 2^j _ 1 V2 . 

4. Computational Aspects 

First of all, it should be emphasized that 
the above analysis is based only on RR them- 
selves without any use of the "integral struc- 
ture" of X^. If there are extra restrictions 
which follow from this structure, they should be 
incorporated into RR. For example, if a loop 
integral admits symmetry under some permu- 
tation of indices, one has to use RR in the 
properly symmetrized form. Thus, in example 
of Sect. 4 X VlV2 = X V2Ul |3J. and if one adds 
RR with permuted indices, then GB becomes 
D\l Vl -i i/ 2 -i) DiD^L Vy -\ ,/ 2 -i, D^Ivi-ivi-i. 

It should be also noted that the GB technique 
is applicable to RR with specified indices. Pro- 
vided with an appropriate ordering on the inte- 
grals (priority criteria |4lfi| ). the GB algorithm is 
just Gaussian elimination Hence, it is con- 
ceptually identical to the Laporta algorithm. 

In the above illustrative and very simple ex- 
ample (Sect. 3) we performed the GB calculation 
over the field. Generally, however, one has to 
be careful in manipulation with arbitrary indices 
occurring in the coefficients of difference polyno- 
mials. In calculating GB over a coefficient field, 
the leading coefficients of difference polynomials 
are treated generically, and division by them can 

tation is very expensive and can be applied only to small 
problems (Sect. 4). 
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be done at the intermediate steps of calculation. 
GB computed this way may loose its GB proper- 
ties under some specifications of parameters. This 
happens when, in the course of reduction, the di- 
vision by a leading coefficient was done, and this 
coefficient vanishes under the specification. To 
avoid this trouble, one can compute GB over a 
coefficient ring rather than over a field (e.g. with- 
out division). But computation over the ring may 
lead to growth of intermediate coefficients and 
make the computation much more tedious. 

As it is well-known (see JT] and references 
therein), the running time and storage space re- 
quired by the GB algorithms the GB tend to be 
exponential and superexponential in number of 
variables (indices of RR) . Besides, the presence of 
parameters and symbolic indices in coefficients in- 
creases substantially the volume of computation. 
Therefore, to be practically applicable to mod- 
ern multiloop calculations, one needs efficient im- 
plementation of the GB algorithms provided by 
special tuning to finite-difference polynomials of 
type ©. 

To our knowledge, there are no GB packages 
available for multivariate finite-difference poly- 
nomials, except those designed for the (noncom- 
mutative) skew-polynomial algebras or Ore alge- 
bras 8 . Our experimentation with the Maple 
package implementing the Buchberger algorithm 
for skew polynomials, revealed its very low effi- 
ciency for RR with parameters. Thus, Maple on a 
2 Ghz PC with 512 Mb RAM was not able to con- 
struct (in reasonable time) GB for the full (but 
still rather small) one-loop example from Z, when 
both m 2 and q 2 are nonzero and present in RR. 
To our opinion, this is because of extra computa- 
tional costs caused by internal noncommutative 
settings of underlying built-in operations. 

We find more perspective to adapt our invo- 
lutive algorithms and software packages [12113) 
to linear multivariate finite-difference polynomi- 
als. The involutive algorithms also compute GB, 
but they are based on the other course of inter- 
mediate computation different from that in the 
Buchberger algorithm. Our implementation of 
involutive algorithms already shown their higher 
efficiency 13 for conventional commutative poly- 
nomials. In addition, involutive algorithms admit 



effective and natural parallelism, that may be cru- 
cial to manage with multiloop calculations needed 
in modern higher energy physics. 
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